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We construct the general 1/2 BPS M2 giant graviton solutions asymptotic to the eleven- 
dimensional maximally supersymmetric plane wave background, based on the recent work of Lin, 
Lunin and Maldacena. The solutions have null singularity and we argue that it is unavoidable to 
have null singularity in the proposed framework, although the solutions are still physically relevant. 
They involve an arbitrary function F(x) which is shown to have a correspondence to the 1/2 BPS 
states of the BMN matrix model. A detailed map between the 1/2 BPS states of both sides is 
worked out. 
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1 Introduction 



Recently, comprehensive works on 1/2 BPS geometries in the type IIB and 11D supergravities [1], 
as well as in the six dimensional supergravity have been done [2-4] . (For other related works please 
see [5].) They are supposed to describe geometries corresponding to 1/2 BPS states of dual CFT 
in the framework of AdS/CFT correspondence [6,7]. The story in the case of type IIB turns out 
to be much interesting, because in this case we know much information on both the geometry 
side [1,8] and the corresponding CFT side [9]. The two dimensional phase space of the emergent 
fermion picture in the 1 /2 BPS sector of J\f = 4 SYM theory is geometrically realized as a surface 
of specific boundary conditions. Detailed matching between 1/2 BPS operators in J\f = 4 SYM 
and the geometries of giant gravitons, in a way suggested by [9,10], provides another convincing 
example of AdS/CFT correspondence. It was pointed out that it is even possible to identify the 
fermion system directly in the supergravity side [12]. 

Contrary to 1/2 BPS geometries of AdS^ x S 5 in type IIB, where we only need to solve a 
linear Laplace equation with a specified boundary condition, the 1/2 BPS geometries in the M- 
theory are much harder to find because we have to solve a continuum version of nonlinear Toda- 
like equation [1]. The fact that we again have a boundary plane divided into regions of two 
different kinds of boundary conditions seems to be similar to the type IIB case, but the nature and 
even the existence of smooth solutions are not very clear yet. However, for the geometries in the 
background 11-dimensional plane- wave, in which case there is a translational symmetry along one 
direction on the boundary plane, it is possible to map the non-linear equation to the linear Laplace 
equation [1,13]. In this work, we analyze this linear equation carefully to find solutions of physical 
relevance. 

The solutions we provide in this work depend on an arbitrary function F(x), which we interpret 
as the number density of 1/2 BPS spherical M2-branes of radius x in the transverse R 3 space of 
the 11D plane wave background. This identification will be obtained from comparing the charges 
of the geometries with those of 1/2 BPS solutions in the BMN matrix model. However, we find 
that the geometries of our interest have null conical singularity. In fact, we will try to give a 
convincing argument that this type of singularity is inevitable in the proposed ansatz of Lin, Lunin 
and Maldacena (LLM). We leave it as an open question to discuss the acceptance of this type of 
singularity. Except this singularity, the resulting geometry seems to have a nice correspondence to 
the BMN matrix model expectations [7] . 

After a brief review on the 1/2 BPS geometries in M-theory in the next section, we construct 
the explicit solutions of 1/2 BPS geometries in the plane-wave background in section 3, followed 
by explicit examples and analysis in section 4. In section 5, we discuss the singularity and argue 
that it seems impossible to avoid the singularity in the given ansatz. We compute the charges of 
the geometries and finally compare them with the BMN matrix model in section 6 and 7. 



2 A Review of LLM Construction in M-theory 

In this section we review the construction of 1/2 BPS geometries in M-theory following Lin, Lunin 
and Maldacena (LLM) [1]. M-theory or the eleven dimensional supergravity in the low energy 
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limit contains two familiar 1/2 BPS objects, M2 and M5 branes. M2 (M5) brane is electrically 
(magnetically) charged with respect to the three from gauge potential of the eleven dimensional 
supergravity. The corresponding supergravity solutions interpolate between flat Mu and AdS^xS 7 
{AdS-j x S 4 ) geometries for M2 (M5) branes, which are known to be maximally supersymmetric 
solutions of M-theory. In addition to this, one additional background is known which is also 
maximally supersymmetric. This corresponds to the plane wave solution obtained from AdS& x S 7 
or AdSj x S 4 by taking the Penrose limit [14]. 

Motivated by AdS/CFT correspondence [6] and its extension to plane- wave/CFT correspondence 
[7], one is interested in the BPS states of the dual CFT obeying A = J, where A is the conformal 
dimension of the CFT fields and J corresponds to the U (1) charge in the R symmetry group. Let us 
consider 1/2 BPS geometries in AdS7 x S 4 . These are associated to the chiral primaries of the (2,0) 
theory. The chiral primaries of (2,0) theory can be described in terms of two dimensional Young 
diagrams. These states preserve half of the supersymmetry and are invariant under the symmetry 
group, SO(3) x SO(6) x R of the (2,0) theory on R t x S* 5 . In the dual supergravity description, 
one is interested in the half BPS solutions with this symmetry. 

The most general ansatz for the eleven dimensional background consistent with SO (3) x SO (6) 
symmetry can be written as [1] 

ds 2 = e 2X (^dnl + e 2A dn 2 2 + ds\ 

G (4) = G pi 1x211311 a dx^ 1 A dx^ 2 A dx^ A dx^ + F MlM2 dx^ A dx^ 2 A d 2 ti 2 2 (1) 

where Q, 2 . and Q 2 being the line elements of the unit two sphere and five sphere respectively and 
being the field strength of the three form gauge potential of eleven dimensional supergravity. The 
Greek indices run over 1, ..,4 and refer to the four dimensional metric ds 2 . A and A are functions 
of all the coordinates and m? is some fixed parameter independent of the coordinates. 

The number of supersymmetries preserved by a background is in one to one correspondence with 
the number of covariantly constant Killing spinors admitted by the background. To find solutions 
which preserve some fraction of the supersymmetry one has to solve the Killing spinor equation [15], 

v m e + (rr r - ^ pqr ) G Wnpq rt = o. (2) 

One can solve the above equation by decomposing the spinor, £ into spinors of ds 2 and seven 
dimensional spinors and perform the dimensional reduction on S 5 . Further reduction on S 2 gives 
the equations for the 4-dimensional spinor. One can fix the form of the metric by working out the 
properties of the Killing spinor bilinears and using Fierz identities. Defining a new coordinate, y 
in terms of the bilinears, the metric components can be written in terms of a single function, D of 
spatial coordinates of the four dimensional metric ds 2 (for the detailed discussion of the reduction 
procedure and fixing the form of metric, we refer to the appendix of Ref. [1]) 

ds 2 = -4e 2A (l + yV 6A ) (dt + Vidx 1 ) 2 + - + C 2 ^_ 6A [dy 2 + e D (dx\ + dxfj 
+4e 2X dn 2 5 + y 2 e~ 4X dn 2 2 
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-6A _ 9yD 



y(i-yd y D) 



Vi = \e l3 d 3 D , (3) 

where i,j = 1, 2. We have not written the expression for the four form Gm, as we shall not need it. 
It is also determined completely from the knowledge of the function D(y, xi,xz)- The parameter m 
in the equation (1) has been set to m = \. The function D(y, xi, x-i) satisfies the three dimensional 
Toda equation 

(d\ + dl) D + d 2 y e D = 0. (4) 

It can be seen from the metric (3) that the coordinate y is related to the radii of the two sphere, 
i?2 and of the five sphere, i?5 by the relation y = i?2-R|/4. At y = either two-sphere or five-sphere 
shrinks to zero size. The boundary conditions at y = are required to be such that the geometry 
remains non singular when either of the spheres shrinks. In the case where two-sphere shrinks while 
the radius of five-sphere remains constant, one finds d y D = at y = and D is independent of 
y. In the other case where five-sphere shrinks to zero size while the radius of two-sphere remains 
constant, we have D ~ logy at y = 0. Thus, one can write the boundary conditions at y = as 

d y D = 0, D = finite, when S 2 shrinks 

y, when S 5 shrinks (5) 



e D 



We shall refer the above boundary conditions as 'finite type' and 'linear type' respectively as the 
function e D becomes finite or goes linearly with y. We shall be interested in solutions of the Toda 
equation (4) satisfying the above boundary conditions. In general the explicit solutions to Toda 
equations are not known, but in the particular case when there is an extra spatial isometry (say x\ 
in our case) the Toda equation (4) reduces to 

d\D + d 2 y e D = 0. (6) 

By a change of variable* 

e D = p 2 , y = P d p V, x 2 = -d v V (7) 
the Toda equation (6) can be mapped to the Laplace equation in three dimensions [13] 

-d P (pd p V) + d 2 V = 0. (8) 

Therefore, the task of solving the Toda equation (6) is reduced to solving the Laplace equation in 
the new coordinates subject to the boundary conditions (5). 



3 M-Theory Giant Gravitons 

In this section we are interested in the solutions of the Laplace equation (8) which correspond to 
giant gravitons in M-theory. (See Ref. [11] for a linearized solution for M-theory giant graviton.) 

*Our notation differs from Ref. [1] by X2 — » —xi- 
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In particular we shall be interested in the solutions which are obtained from superposition of the 
solutions of the Laplace equation. The plane wave background is given by a solution of the Laplace 
equation [1]; 

V = p 2 V~lv 3 ■ (9) 

As a warm-up exercise, let us discuss the boundary condition for this solution explicitly The metric 
for this solution can be written as 

ds 2 = -A{p 2 + Ai 1 2 )dt 2 -Adxidt + A(dp 2 + dif^ + Ap 2 dtt 2 5 + Ar] 2 dQ 2 . . (10) 

We see that p and r\ are nothing but the radii of the five-sphere and two-sphere respectively. We 
are interested in the issue of boundary conditions at y = in terms of the new coordinates (p, rj). 
They are related to the previous (y, x-i) coordinates by 

y = 2p 2 v , X2 = _( p 2_ 2r?2) _ (n) 

The boundary at y = in the (p, ry)-plane corresponds to either p — ► (S* 5 shrinking), while keeping 
rj (radius of two sphere) finite or rj —> (S 2 shrinking), while keeping p (radius of five sphere) finite. 
When the five sphere shrinks one can see that the boundary conditions at y = are of the 'linear 
type' that is e D ~ y. In the other case when the two sphere shrinks, the boundary conditions at 
y = are of the 'finite type' that is e D ~ finite. 

Now we shall construct solutions which are superpositions of the solutions of Laplace equation 
(8) upon the plane wave background. Later we shall identify these solutions as half BPS states of 
BMN matrix model [7]. Solutions to the Laplace equation have several types. There are polynomial 
type solutions like the case of the above plane wave solution. To recover the plane-wave background 
asymptotically they should be absent. Other class of solutions are the superposition of 

K (kp)e lk \ I (kp)e ik \ (12) 



or 



Hkp)e k \ N (kp)e k \ (13) 

in terms of the Bessel functions. After careful analysis of these types of solutions, it turns out 
that the relevant one for our application is the solutions involving Ko{kp), the zeroth order Bessel 
function of the imaginary argument, which shows the logarithmic singularity at p = 0. For this 
choice, generic solutions to the Laplace equation (8) can be written as 

Vi = 2 f°° dkA{k)K {kp)e ikr! . (14) 

Using the integral representation for the Bessel function Kq(x) 

/ x f°° , costa 1 [°° , e itx 

the final form of the solution obtained from superposing on the plane-wave background is 



V = V + Vi = p 2 7] - -rf + 



A 



where the function F[x) is the inverse Fourier transform of A(k). In fact, V\ is the electromagnetic 
potential of a line charge density F{rj) distributed along the ry-axis. 



The other choices of polynomials does not lead to the asymptotically plane wave geometries. 
One may also show that the choices of Iq, Jq and iVo are not consistent with the asymptotically 
plane wave geometries with the desired boundary conditions. In section 5, we will actually consider 
the most general acceptable form of V\ for the solutions of the Laplace equation, and compare with 
the properties of the solutions in this section. 

Using (7), the expression for y and xi may be written as 

■2 



OO -j- 

y = 2p' V + p I dt-=— F'(r] + tp) 

oo V 1 + t z 



x 2 = -p 2 + 2?7 2 - 



50 , 1 

dt . 

oo y/1 + t 2 



F'(V + tp) 



(17) 



where prime denotes taking derivative with respect to the argument. 

In terms of the coordinates (p, rj), the expression for metric can be written explicitly as: 

K 



ds< 



o6A 



-4A 



4e 

+Ae^dQ'i + y l e 



W (dt + V^f + ^e 6X (d Xl ) 



+ 



= W- 



Vl 



K drj 



(18) 



where W 



and K = pj, J being the Jacobian of transformation from (y, x-i) coordinates 
to (p, rj) coordinates (see below). Unlike the case in terms of the coordinates (y, X2), the solution 
here is no longer implicit. Namely, all the metric components as well as the field strengths can be 
written explicitly in terms of F(x), though the actual expression may not be much informative . 

Let us discuss the boundary conditions for these solutions carefully. To get the right kind of 
boundary conditions (5) for giant gravitons, we restrict the function F(x) to be an odd function. 
This choice will be further explained in the later section. In addition, F(x) should be nonnegative 
for x > 0. (For the negative F(x), the y = line stops somewhere in the r\ axis. In this case, at the 
boundary of r\ axis excluded by the y = line, the radii of two and five spheres become finite while 
p goes to zero. For a smooth geometry, one has to utilize x\ circle but there is no way to match 
the period of the x\ circle except for some special cases that are not of our interest.) Other than 
these restrictions, F(x) can be quite general; F(x) may even be discontinuous as will be illustrated 
in the next section. 



Now let us consider the union / of intervals, /„ = [a n , b n ] with < a n < b n , which is ordered 
and non overlapping. Namely b n < a n+ \. If there are I such intervals, / is given by 1\ U I2 U • • • U J/. 
Without loss of generality, one may consider F(x) being nonzero only in x € /. For this setting of 
F(x), the boundary at y = takes in general the following shapes in the (p, ry)-plane. First it is 
clear from (17) that 77 = solves y = as F(x) is chosen to be odd. Furthermore, p = solves 
y = equation if r\ does not belong to /. Thus one may see that p = 0-line and r] = 0-line together 
form some part of the relevant y = boundary. The coordinate (p, rj) is defined in the region 
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within the first quadrant, i.e. p > and rj > 0. Around each interval I n , there is a larger interval 
I n (I n C / n ), in which y = boundary is given by p = g n {v) > 0- The intervals /„ may overlap 
with each other generically. If this happens, the overlapped intervals are merged and we may get 
a new set of interval {I n }- This new set is constructed such that each interval I n is disjoint with 
each other and that p = g n {v) / within each interval. It is also possible that the first interval 
may be extended to the p axis, which means p = g\(0) > 0^. In conclusion, the y = boundary 
consists of p axis with p > gi(a\), rj axis that does not belong to 7, and the collection of p = g n {rj) 
lines for rj £ I. We illustrate the shape of boundary for I = 3 case in Fig. la by bold line. 



) 

) 
) 

I P 

a) y = boundary b) y = boundary in (1, 2) plane 

Figure 1: An illustration of the shape of y — boundary in (p, ?y)-plane and in (x\, x 2 ) for I = 3. 

Along the boundary in (p, rj) space, that is, along the y = level curve, X2 has a monotonic 
behavior. For instance, X2 in Fig. la is monotonically increasing along the boundary in the direction 
of the arrows. The proof of this follows from the fact that the Jacobian is positive definite for y > 0, 
which will be shown shortly, and that y is increasing to the right of the boundary. Then for each 
interval I n , one has corresponding intervals I n = [a n , f3 n ] in the X2 coordinate. From the monotonic 
behavior, the intervals T n 's are again ordered and not overlapping with each other. In these intervals 
as well as the X p i ane = [-co, with (3q = X2(rj = 0, p = g\(ai)), the finite type boundary condition 
is realized. X p i ane may be viewed as M2 branes responsible for the plane wave background. The 
intervals I n in the X2 direction is considered as representing extra M2 giant gravitons. The shape 
of boundary in the [x\, X2) plane is depicted in Fig. lb for the I = 3 case. The number of M2 giant 
is not given by the area of the strip but given by the area weighted by p 2 as will be discussed in 
the section 6. 




Let us give more details for I = 1 case where F(x) is non zero only in an interval [a, b]. Integrating 
by parts and Taylor expanding the function F(x), we can write the expression for y (for small p) 
as: 



where 



y = -2F( V ) + P - [d v x 2 (p = 0) + 3p 2 g( V ) + .. 
9(V) = f 

J —1 



(19) 



roo F(rj + w) 
aw — : — 



fNote that ai is zero in this case. 



6 



Prom this expression we note the following: 

i. For p = and rj £ [a, b], one can see that y ^ and hence p = and r/ £ [a, b] is excluded 
from the boundary. 

ii. The condition that y > for infinitesimal p further requires x' 2 > and one can see that the 
intervals [ry m j n ,a] and [b,r] max ] are also excluded from the boundary, 7? m i„ and 7? maa ; being the roots 
of d v X2 = to the left of a and right of b respectively. 

iii. The p > and r] ^ [rimin,r]max] region is also not a boundary except for r\ = 0, which is a 
trivial boundary with boundary conditions of finite type 

p 2 = e D ~ finite 



iv. For p = and 77 ^ hmim'/mm] the boundary conditions are of linear type as the function 
F{rj) vanishes outside the interval [rjmimVn 



\max\ 



2 D 
p = e ~ y 



v. For p > and rj G [r] m i n , i] max ], the boundary conditions are of finite type as y = boundary 
looks like a bump toward p > region 

p 2 = e D ^ finite 

However, we always find that <9 y .D = condition is not satisfied in this part, and the geometry has 
a null singularity as will be discussed in section 5. Thus a suitable choice of the function F{rf) (odd 
and positive in a finite interval) gives the solutions which obey the almost right kind of boundary 
conditions for a giant graviton except the appearance of null singularity. In the next section we 
shall discuss this in an explicit example. 

By now we have not worried about smoothness of the coordinate transformation from (y,X2) 
coordinates to (p, rj) coordinates. For a well behaving solution in (y, X2) coordinates one would 
require that the change of coordinates should be well-defined throughout the region. To analyze 
this, consider the Jacobian of the transformation (y,^) — > (p, rj) 

j = dy_ 8x2 _ dy 8x2 , . 

dp drj dri dp ' 



Using the definitions (7) of y and X2 and the equation (8) one finds 

1 



J = 
P 



dy\ 2 (dy 



dp) \drj 



(21) 



which is positive semi-definite. 



Also we have implicitly assumed that y > 0. Since y is related to the radii of the two-sphere and 
five-sphere, one would like to show that radii of two-sphere and five-sphere never become negative. 
The radii of two-sphere and five-sphere are related to the quantity e 6A defined in the metric (18). 
The requirement for the non-negativity of e 6A reads from (18) as 

K > 2-$ (22) 
pdp 
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which is true if the following inequality holds 



dy 2y 
dp - p' 

Integrating by parts the expression (17) for y becomes 



/CO 
dx 
-co 



Fix) 



(p 2 + (x - n) 2 ) 2 



(23) 



(24) 



Taking partial derivative with respect to p and using the fact that F[x) is an odd function, one 
gets 



dy = 2y + 
dp p 



3p / 

J o 



dx 



P 



F{x). 



(25) 



ip 2 + (x — 7]) 2 ) 'i ip 2 + (x + 7/) 2 )2 " 

^From this, the inequality (23) follows provided F(x) > for x > 0. Hence, the non-negativity of 
y further requires that F(x) > for x > 0. 



4 Examples of Giant Gravitons 



To give a compelling basis to the validity of the solutions we have described, let us analyze an 
explicit example more carefully. For this purpose, we make the simplest choice of the function 
F(x), that is, the one with <5-function distribution, 

Fix) = f (6(x - no) - S(x + no)) . (26) 

It is easy to calculate y and X2 to find 

y = p 2 (zv-fo— — — ^ — r^r + /o— — — ^ — — t) , 
[p + {v - vo) ] 2 [p + yn + vo) } 2 

x 2 = -p +2r/ +/ — — ^TT-/o— — ; — s • (27) 

[p z + in - no) ] 2 [p + iv + vo) \ 2 

We are interested in the region of p, r/ > in which y > 0, and the y = boundary is of special 
importance as discussed before. One obvious component of y = curve in (/?, r/) plane is r/-axis 
{p = 0), whereas there may be other components by solving 



2r/ = /o( 



" ^Tt) • (28) 



[p 2 + in~ iio) 2 ] 2 [p 2 + (11 + vo) 2 } 2 



We will find in the following that there are indeed two additional components from this equation; 
one is the p > 0-axis (n = 0) and the other is an approximate semi-circle (or "bump") around the 
point ip, n) = (0, 770) in the p > region. Let us call the two intersection points of this approximate 
semi-circle with the ?7-axis, A and B, whose positions are (0,r/ m i n ) and (0, rj max ). At first sight, this 
may seem to indicate a problem because y = 0-curve bifurcates at A and B into p = 0-curve and 
the semi-circle. However, we will show that, for n £ [f]mm,Vma,x], we should actually discard p = 
segment and take only semi-circle part. This is because the semi-circle turns out to be the correct 
boundary of y > region, while p = 0, 77 £ [rjmin, ijmax] is the boundary of y < region instead. 
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Therefore, the whole picture of y = boundary for y > region would be 

{p > 0, rj = 0} U {p = 0, 7] e [0, r] min ] U [r/max, oo]} U {p = #(77), 77 G [77 min , r/ max ]} , (29) 

where in the last component, p = g{rj) is the equation of the (approximate) semi-circle, whose 
implicit expression is given by (28). In the previous section, we have shown that Jacobian between 
(p, rj) and (y, x-i) is positive definite. This implies that %i is monotonic along the level contour 
of y = in (p,rj) plane. Recalling the map e D = p 2 , we thus have two separate components of 
D ~ (finite) boundary condition as y — ► 0, that is, the first and the last component in (29), while 
the second component in (29) corresponds to the boundary condition e D ~ cy as y — ► 0. Therefore, 
this solution is describing one-band excitation of the plane wave background, whose interpretation 
would be a supergravity realization of 1/ 2 BPS M2-brane. 

We can explicitly confirm the claimed features of y = boundary of y > region in (p, rj) 
plane by considering sufficiently large 770 S> 1 with /o ~ 0(1). We may reasonably expect these 
properties (e.g. the existence of the "bump") not to change for generic 770 and /o, indicating a large 
degeneracy of solutions for a given topology of one-band excitation. In solving (28), the solution 
rj = is obvious, and we look for solutions with 77 7^ 0. For large 770, the right hand side of (28) 
will be small unless 77 ~ 770, and the left hand side tells us that rj is small, allowing us to expand 
the right hand side near 77 = 0. This in fact gives us back the 77 = solution. However, near 77 ~ 770 
region, we may keep the first term in the right hand side while the second term is neglected. This 
gives us the equation 



2 2 
3 



describing a semi-circle of radius {J^j^ centered at (0,770), and the constants 7? m i n / max is given by 

»±(A)'- 

Around this semi-circle, we have 

■ (3i) 

from which it is clear that the semi-circle bounds the y > region to its right, while p = 0, 
V ^ b/mim f?max] segment is the boundary of a wrong region. It is also seen that for rj £ [0,77 m i n ] U 
[r/ max , 00], the p = line should be taken as the boundary of y > as in (29). What remains to be 
shown is that p > 0, 77 = is the boundary of the y > region to its above. Near 77 = 0, we have 

y*2p*Jl- 3for >° , )+0(pW) ■ (32) 

Therefore, by taking 770 and /o such that (?7o) 4 > 3/o, we have y > above the rj = boundary. 
Moreover, the arguments in the previous section guarantee that y > remains true in the whole 
region bounded by (29) and there is no other boundary. We hope it is by now convincing that 
(29) is the correct boundary of our solution, and it represents an example of 1/2 BPS one-band 
excitation of the plane wave background. 

It is straightforward to check that these solutions have the correct boundary conditions as we 
approach y — > 0, that is, e D = p 2 behaves either ~ cy or (finite), except the segment {p = g(i]),rj G 
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[ijmin, Vmax] } on which d y D / 0. Therefore, we have a band of conical singularity which turns out 
to be null-like. In seeing this, it is useful to write the resulting metric explicitly after we change 
the coordinate from (y, x<i) to (p, 77) as in (18). The reason for this is because the presumably 
complicated function e D of (y, x-i) is simply p 2 in (p, rj) plane. Moreover, (p, 77) is preferred for 
another reason; they become the usual radial coordinates of R 6 and R 3 in the asymptotic plane- 
wave background, and this fact will be useful when we calculate charges of our solutions in the 
later section. 

It is easy to generalize the above discussion to find explicit multi-band solutions. By simple 
superposition of one-band solutions, we have 

y = 2 V p 2 + P 2 j2{-fo r - - 1 ^r + /o r : ; 1 > 

1 [p 2 + fa - »#) 2 ] * [p 2 + (r] + »7o) ] 5 

X2 = V + 2r ? 2 + E{ + /o ^ r~fS ^ A ■ (33) 

For sufficiently far separated t]q and small ftf, each n-th term independently generates small semi- 
circle bump around 77 ~ tjq in the y = boundary of y > region. The component of {p > 0, 77 = 0} 
and {p = 0, rj € i? + — U^it^mim ^max]} are as before. Because the behavior near each bump is 
not affected in an essential way by other bumps, the properties of these solutions, like the nature 
of singularities, should be same as in the one-band solutions, at least for generic values of t]q and 
fo". As it will become clear in the later section, the charges for multi-band excitations turn out to 
be additive, conforming to the naive expectation. 

We would like to mention that the above solutions obtained from 5 function distribution of F(x) 
are mere examples of vast possible choices of F{x). For instance, we may try a step function: 
F(x) = fo for x £ [a, b], —fo for x G [—6, —a], and otherwise zero. We have analyzed this case to 
find one-band excitations for reasonable values of fo, a and b, and also expect multi-band solutions 
by superposition of them. We do not believe that different choices of F(x) are related to each other 
by coordinate reparametrization, and this suggests a huge degeneracy of solutions parameterized 
by a single function F(x). The physical interpretation of this function will become clear in section 
6 when we discuss charges of the solutions. 



5 Null Singularity 



Our solutions presented in (17) are singular at y = if i] G 7. In this region of y = with rj E I, 
the finiten 
Note that 



the finiteness condition of D is satisfied but d y D = is not satisfied. This can be seen as follows. 



d y p 2 ~2d p y V 

where the first equality is from (3). From the second equality, one may see that d y D = 2d y hip = 
condition is equivalent to d p y = 0. However one may see from (25) that d p y is finite (and positive) 
for a finite p > and for any nontrivial F (F(x) > when x > 0). Thus the boundary condition is 
not satisfied for 77 6 I leading to the above mentioned singularities. At the singularities, one finds 
that e 2A as well as the radii of both the two and five spheres are vanishing. These are null-like 
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conical singularities and not much different from those of extremal branes, whose horizon is singular 
in general. 



Before we discuss the nature of the singularities further, we like to argue that these singularities 
are in fact unavoidable in our framework for the geometries of 1/2 BPS giant gravitons asymptotic 
to the plane wave geometry. The most general solution of the Laplace equation that is asymptotic 
to the plane wave geometry is given by 



x> roo r2n 

dx / dss / d9 

oo JO Jo 



F(s,x) 



\J p 2 + s 2 — 2s p cos 9 + (rj — x) 2 



(35) 



where the source F(p, rj) should be distributed outside of y > region in (p,r])-p\aiae. In writing 
the last part V\ in (35), we used the fact that V\ satisfies the three dimensional Laplace equation 
with the rotational symmetry around rj axis, and also the boundary conditions at infinity for the 
asymptotically plane wave geometries. Let p s be the maximum radius of the nonvanishing source 
distribution occurring at rj = ri s . After a short calculation, we find 



Op P J-oo Jo Jo 



-2s cos 9F(s, x) 3(p 2 + s 2 - 2sp cos 9)F(s, x) 



_3 
P Z2 



+ 



, (36) 



where Z 



p 2 + s 2 



2sp cos 9 + (ry 



As before, requiring F(s, x) be odd in x and positive 



definite for x > for the right shape of the y = boundary as in Fig. la, it can be easily proved 
that the second term in the above is always positive. This implies that y remains positive in the 
bulk, and we have conical singularities on the "bump" of y = 0-boundary, exactly same as in the 
line charge case. 

Indeed, the potential (35) is equivalent to the one from a line charge density for outside observers 
with p > p s . Using the definition of Kq, the V\ part may be rewritten as 

V 1 = — [°° dk 1°° dx 1°° dss e^^F^x) d9K (kZ(9)) (37) 
2ir J-oo J-oo Jo Jo 



where Z 2 {9) = s 2 + p 2 — 2spcos9. We then perform the 9 integration. Using the formula 

1 C 27T 



1 f Z7r 

— / d9K {kZ{9)) = I {ks)K (kp) for s < p , 

2ir Jo 



one may rewrite (37) as 



Ui 



/_ 



Fix) 



oo y/p 2 + {rj - x) 2 



for p s < p, 



where we define 



/oo roo roo 

dk dy dsse ik(x ~ v) I (ks)F(s,y). 
-oo J —oo Jo 



(38) 



(39) 



(40) 



Therefore, for the charge distribution localized around the r\ axis with the cylindrical symmetry, 
the system may be equivalently described as a line charge density at r\ axis for the observer outside 
the maximum of p (p s ) of the charge distribution. If F(s, x) is odd under x, F(x) is also odd. The 
nice properties derived in the previous paragraph can be seen more easily in this context because 
we showed in section 2 that they hold for line charge distributions. 
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We think there is in fact a stronger argument for the singularity on the "bump" of y = 
boundary. Quite generally, suppose we have a solution with the necessary boundary condition on 
the y = curve satisfied, and imagine focusing on the edge point where the boundary condition 
changes from the linear type to the finite type. Note that for this purpose, it will be sufficient to 
assume that we can go to the (p, rj) description only locally around this point, and the conclusion 
we draw from this is universal in this sense. In (p, r/)-plane, the y = curve leaves the 77-axis at 
this point, and let us consider the expansion of the potential around it. (This kind of point should 
exist for the giant gravitons unless y = line stops at the rj axis.) The general argument shows 
that this kind of point should occur away from the charge distribution. 

Taylor expansion around the point (0, 770) should give us polynomial type of solutions of the 
Laplace equation, which can be constructed order by order by the formula 

V = Y^Fn( P 2 W , (41) 

n>0 

with the recursion relation 

Fn + 2(x) = ~ {n + 2) 4 (n + 1} 0<(*) + KW>) > ( 42 ) 

where fj = 77 — t]q. The sum is actually finite for each independent solution. To describe the y = 
curve leaving 77-axis, the lowest order expansion of the potential V looks like 

2 8 
V = co + afj + c 2 (p 2 7? - -f) + c 3 (p 4 - 8pV + -t) 4 ) + • • • (43) 

Then y = pd p V becomes 

y = 2p 2 (c 2 f/ + 2c3p 2 -8c 3 rf + •••) (44) 

where c 2 should be nonzero for the boundary condition of linear type. One may see that y = 
curve is described by p = and the curve of 

c 2 ?? = -2c 3 p 2 + • • • (45) 

The evaluation of d p y leads to 

d p y = Ap(c 2 fj + 4c 3 p 2 - 8c 3 f/ 2 • • •) (46) 

Thus along y = with p / 0, we see that d p y / in general. In case c 3 = 0, one has to consider 
the higher order terms but the conclusion is not changing, unless V(p, rj) exactly factorizes into 
F(p 2 )G(r]) where we find no relevant solution. 

As mentioned earlier, the only remaining possibility is that y = line stops at somewhere in 
the 77 axis. This is possible as y works as a coordinate. A possible example arises for the case of 
F(x) < for x > 0. In this case, at the boundary of 77 axis excluded by the y = line, the radii 
of two and five spheres become finite while p goes to zero. To make a smooth geometry, one has 
to utilize x\ circle but there is no way to match the period of the x\ circle except for some special 
cases that are not of our interest. 

This completes our argument for the non-existence of regular 1/2 BPS giant solutions that are 
asymptotic to the plane wave geometries. Therefore in the supergravity description of the ansatz 
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we are considering, the occurrence of singularities seems inevitable. There are many examples of 
singular solutions in supergravity theories. As mentioned earlier, the singularities of our solutions 
are null-like and rather similar to the null singularities occurring in the supergravity description 
of extremal D-branes. In case of extremal D-branes, we view the D-branes as extra fundamental 
degrees of freedom corresponding to sources of supergravity fields. Without the source terms 
describing the D-brane dynamics, the supergravity description would miss some part of the degrees 
of freedom leading to singularities. Indeed the structure of singularity matches precisely with 
the Born Infeld type source term for D-branes. Whether M2-brane description [18] of the giant 
graviton dynamics matches with the singularity structure of the above solution is of particular 
interest. Further studies are required in this regard. 



6 Charges 

In this section, we like to identify the charge content of the 1/2 BPS system. This will be essential 
for the later comparison with the 1/2 BPS states of the matrix model in the later section. First 
consider the lightcone momentum, p + the charge associated with the translation along x\ direction, 
which counts the number of DO branes in the IIA string theory. 

Since the solutions are asymptotic to the plane wave, the usual ADM definition of energy and 
momentum will not work for our case. In Ref. [16], a method of finding the conserved quantities 
for the background with translational isometry was discussed. Taking x\ as a relevant isometry 
direction, one may get the expression for the charge straightforwardly. Here we would like to present 
an alternative method for computing the charge, which gives precisely the same results. We first 
go to the IIA string description compactifying the eleven dimensional theory on the circle along 
x\ direction. The ten dimensional metric in the string frame is related to the eleven dimensional 
metric by 

ds 2 n = e~¥dsj + ef {dx\ + C^f (47) 

where 4> is the dilaton field, is the coordinate for the ten dimensions and C M is identified with 
the R-R one form potential. The integral of *dC over the eight sphere at infinity corresponds to 
the charge. Working with the following general form of the metric, 

ds 2 = -h t dt 2 + 2w dtdxi + h x dx\ + h ij dy t dy j , (48) 

the IIA metric becomes 

ds 2 10 = ^/K dt 2 + hijdtfdy^ , (49) 

where we introduce L = w + hfh x , and e 3 = h x . The only non- vanishing component of C M is 
Ct = w/h x - The expression for the charge reads 

Q = _L [ *dC = — g- [ Vdeth [-dih x + h x d iW /w] (1 + h t h x /w 2 )- 1/2 h ij nj r 8 dn (50) 
11 11 Jy = oc 

where r 2 = y l y l (~ p 2 + rj 2 for large r), the integration is over unit eight sphere and we restore 
2kfi dependence. The vector m that is normal to the eight sphere is defined as nini = 1. 
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Let us consider the solutions with general F(x), for which the metric components at large r 
behave as, 



ht = 0(r 2 ) , w = -2 + 0(l/r 5 ) , h l3 = <% + 0(l/r 5 ) , 



and 



h 3 



2yp 



pdy 



ydp 



3p jjg 

r/r r' 



where p and q are defined by 

^ /"oo J roo 

p = - / dxF(x) , q = - dxxF(x) . 

^ J — oo ^ «/ — oo 



(51) 
(52) 

(53) 



The first relation in (52) follows directly from the metric (18). The asymptotic expansion of h x 
may be obtained from the expansion of y in (24), 



2p 6q 



(54) 



From the asymptotic behavior of the metric components, it is clear that Q diverges once p is 
non-vanishing. This is one reason why we choose the function F(x) to be odd. We do not claim 
that the solution is unphysical but only that the charge p + is not well-defined for nonzero p. For 
our choice of odd function F(x), the charge Q becomes 

105 uj 8 

where denotes the volume of the eight sphere of unit size. 

Now let us discuss the number of the M2 giants. As discussed in Ref. [1], this may be obtained 
by evaluating the flux *Gu\ through seven cycle (£7) of the five sphere fibered over two-surface 
which ends on the y = region where the five sphere shrinks. This leads to 

W „ f , 



N 2 / *G (4) = / dx 1 dx 2 2e u \ y=0 , (56) 

JY.7 JV 

where V is the region in the 12 plane where S 2 shrinks. The evaluation of this number for general 
F(x) appears to be complicated. But when F(x) is small enough and non-vanishing only for large 
x, N2 is given simply by 

roo 

AN 2 ~ / dxF{x) , (57) 





where A./V2 counts only the extra M2 branes. When F is small enough, the main contributions 
comes from the intervals /, for which 2e D ~ F(j])/r) and dx 2 ~ 4rjdr). Note that the argument x in 
F(x) is related to r\ as one may see for instance in the expression (24). Since 77 is closely related to 
the radius of the two-sphere, the natural interpretation of the variable x would be the size of giant 
M2-brane. In comparison with the BMN matrix model in the next section, this interpretation will 
be clearer. Since the integral Jq°° F{x) dx counts the number of M2-giants, F{x) may naturally 
be interpreted as the number density of the giant M2 with respect to its radius x. Since p + is 
proportional to J °° F(x)xdx, one can see that each M2 should carry the p + charge proportional to 
the size x. Later in the BMN matrix model, one may indeed confirm this prediction. 
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7 Correspondence to the BMN Matrix Model 



The dynamics of the plane wave geometry has a matrix model description. The matrix model may 
be obtained by the standard method of regularizing the supermembrane action in the maximally 
super symmetric plane wave background [7,17]. The x~ (x\ in the previous section) is the isometry 
direction and one is compactifying this direction on a circle of radius R and x + (t in the previous 
section) serves as the lightcone time direction. The conserved momentum of the supermembrane 
along the x~ direction is quantized p + = p- = N/R. In the IIA description, the integer N counts 
the number of DO branes. The matrix variables are given by N x TV Hermitian matrices as they 
represent the strings connecting N DO branes. 

The BPS states of this BMN matrix model have been studied in detail [7, 17, 18]. The 1/2 
BPS states corresponding to the giant gravitons are particularly simple. The bosonic part has nine 
spatial directions described by N x N hermitian matrices and has 50(3) x 50(6) R-symmetry. 

The 1/2 BPS states are governed by 

[X a ,X b ] = ie abc X c (58) 
where the indices a, b, c run over 1, 2, 3 and the other directions remain unexcited. 

The construction of the general 1/2 BPS states is simple. Let us first take an n dimensional 
irreducible representation of the SU(2) algebra where n is a positive integer. The Casimir I? = 
X 2 + Y 2 + Z 2 is given by j(j + 1) = (n 2 - l)/4 where j = (n — l)/2. Thus the radius of the fuzzy 
sphere corresponding to the giant graviton is given by 

R, = ^- <») 

Then generic 1/2 BPS states are given by any combination of the irreducible representations, 

T = [m] © [n 2 ] © [n 3 ] + • • • © [n k ] (60) 

where J2i=i n i = N and we take n\ < n 2 < ■ ■ ■ < n^. This may be thought of as representations of 
two dimensional Young diagrams with total N boxes. Since the gauge equivalent representations 
are not distinguished, we only care about the combinations and the ordering of the dimensions does 
not matter. The 1/2 BPS state T may be presented alternatively by 

oo 

T=J2®fn[n] (61) 

n=l 

where /„ is the multiplicity of the n dimensional irreducible representations. The sum Y^=\ fn 
counts the number of total giant gravitons. 

The dimension n of the irreducible representation is related to the radius of the fuzzy sphere 
and it is natural to relate it to the argument x of F(x), the function characterizing our solutions. 
One may then see the correspondence 

r 'dxxF(x) fn/„, (62) 
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which counts the number N of the DO branes. In addition, one has 



o 



oo ^ 

dxF{x) <- £/ n , (63) 



n=l 



which counts the number of the giant gravitons. Therefore, the physical meaning of F(x) is the 
density of SU(2)-representations of dimension x in large x continuum limit. 

This completes the description of the correspondence. One can see that the data for the giant 
gravitons in the supergravity description is encoded in the arbitrary function F(x) for the interval 
x > 0. The arbitrariness of F(x) corresponds to the arbitrariness of the number of giant graviton 
for a given dimension n. Hence one may say that the correspondence works for any possible states. 
Namely any 1/2 BPS states of the matrix theory specified by {/i,/2,/3 ■ ■ ■} has a corresponding 
supergravity data specified by the function F(x) for x > 0: 

{/i,/2,/s---} <- F(x)forx>0. (64) 

Thus the counting of states for a fixed charge p + is well defined in the matrix model description. 
For the case of the supergravity, the function F(x) is continuous and the volume of the fluctuation 
of F(x) for the fixed p + = J °° dxxF(x) will become zero because we are dealing with infinite 
dimensional space. The situation here is rather similar to the counting of the fluctuation of the 
supertube moduli space [19]. The trouble follows from the continuum nature of the description; 
What we need is a quantization or regularization. 

Such a regularization is already achieved in the matrix model description. The state are even 
protected by the quantum corrections due to the enough number of supersymmetries. The counting 
of states for a fixed p + = N is simple. It is the counting of all possible combination of the states 
{/lj /2, /3, • • •} while fixing the sum J2™=i n fn- This problem may be mapped to the problem of 
the bosonic oscillators with Hamiltonian, 

oo 

H na n a n (65) 

n=l 

where a n satisfies the commutation relation [a n ,ajj = 5 nm while all the remaining commutators 
vanish. Then /„ corresponds to the occupation number of the n-th oscillator a n . The partition 
function is given by 

oo 

Z(w) = tie~ pH = H (1 - u/T 1 (66) 

n=l 

where w = e~@. It is related to the Dedekind eta function 

oo 

V ( T ) = e i7TT / 12 Jl(l- e 2ninT ) (67) 

71=1 

The degeneracy of the state for the fixed N is obtained by 

f Z(w) dw , . 

*->5**S5- (68) 

Evaluating this using the saddle point method, one obtains 

T. 



d N ^V^, (69) 
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for the large N. The entropy is then given by 



S N = lnd N ~ K^j-N. (70) 

The counting problem here basically corresponds to finding the possible ways of partitioning N 
into positive integers. 

Since our starting representation in (60) is the representation of the 2d Young diagram, the 
above chiral boson problem is equivalent to the problem of the 2d Young diagram. Or the above 
(chiral) boson may be mapped into the NS sector of a complex fermion through the well known 
bosonization technique. Then the states of NS sector of the fermion is in one to one correspondence 
with the 2d Young diagram [20]. This 2d Young diagram may be mapped into the c=l matrix 
model with a harmonic oscillator potential, which is related to the N noninteracting fermions 
under a harmonic oscillator potential [9]. But from the data of the supergravity, this relation to 
the fermion seems unmotivated. Namely, we do not find any shape of Fermi liquid in (xi, X2) 
plane of the supergravity description unlike case of the 1/2 BPS geometry of the type IIB string 
theory [1]. 



8 Conclusion 



The giant gravitons in various AdS n x S m spaces as well as in the plane wave limit have been 
studied in the context of AdS/CFT correspondence. Recent progress on constructing full back- 
reacted supergravity geometries of 1/2 BPS giant gravitons, initiated in Ref. [1], gave us much 
insight on the correspondence including the emergent fermionic system in the type IIB string theory 
in AdSs x S 5 . However, given the situation where smooth solutions in the asymptotic plane- wave 
background in 11D supergravity have not been found, it is an important problem to find solutions 
in the recent framework of LLM in Ref. [1]. In this work, we performed a careful analysis for this 
problem, starting from the linear Laplace equation that Ref. [1] proposed. We found a class of 
solutions that involve an arbitrary function F(x) whose interpretation has a nice correspondence to 
the BMN matrix model, that is, it is the number density of spherical M2-branes in the transverse 
-space with respect to its radius x. However, we found that null-like conical singularities are 
present along the M2-branes, and tried to argue that this is a necessary consequence of the given 
ansatz. 

There is one possible scenario to understand the nature of troubles that we encounter in this 
paper. Let us first note that analogous problems arise in the study of charged rotating BPS black 
holes in four dimensional N=2 supergravity theories. All the known solutions have naked singular- 
ities. The dual descriptions by strings or D branes are well understood and counting degeneracy 
is straightforward [21]. Thus the corresponding rotating supergravity black hole solutions should 
exist. In Ref. [22], it is shown that nonsingular BPS solutions may be constructed by the addition 
of massive Kaluza Klein fields in such a way that the solutions decompactify near the core to five 
dimensional black hole solutions with regular horizon. Only asymptotically, they approach four 
dimensional geometries (times circle). If one tries to eliminate the massive Kaluza Klein modes 
by smearing over the circle, the geometries become the known four dimensional solutions that are 
naked singular. If this phenomenon happens in our case too, it implies that the assumption of the 
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isometry along x\ circle is too strict. Namely, dropping the ansatz of the isometry and solving 
the full three dimensional Toda equation, there may be a possibility of resolving the singularities. 
In this case, our solutions clearly have their physical meaning and validity once away from the 
singularities; Only the core singularities need to be remedied by uplifting to higher dimensions. Of 
course our problem is then naturally extended to the problem of finding 1/2 BPS giant gravitons 
asymptotic to AdS± x S 7 or AdS-j x S A geometries. Further studies are required in these directions. 
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